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Abstract 

The radiative decay B* — > B q ^{q = u,dors) is reexamined with a modified light-cone QCD 
sum rule method, in which adequate chiral operators are chosen as the interpolating fields in the 
correlators used for a sum rule estimate of the relevant coupling gB*B q -y- The resulting sum rules 
not only show the physical picture consistent with the underlying physics in B* — > Bgj but also 
avoid the pollution by the nonlocal matrix element (7(5) \q {x) 7^75? (0)| 0), which starts with 
twist-3 and thus may bring a large uncertainty to the sum rule predictions. Also, a comparison is 
made with the previous results from light-cone QCD sum rules and chiral perturbation theory. 
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1. INTRODUCTION 



The physics of exclusive B decays is particularly interesting for a detailed test of the standard 
model (SM) and a careful search for potential signature of new physics. With a large amount of 
available data from the BaBar and Belle, it is believed that a considerable progress will be made 
in the area in the near future. A good many theoretical works contribute to this subject. A recent 
progress deserving mention is the presentation of the QCD factorization formula for B — ► ttti, 
7rK and 7r£)[Q]. Nevertheless, to make a reliable estimate of physical amplitudes demands that 
we calculate precisely the relevant hadronic matrix elements. It is practically beyond our present 
ability to do such a calculation from first principle. At present, most of efforts have been devoted 
to looking for phenomenological approaches to nonperturbative QCD dynamics. Among the most 
successful models is QCD sum rule method [T2J. As it is often case, however, this approach has 
a restricted application to heavy-to-light transitions, failing to give a viable behavior of the form 
factors with the heavy quark mass at small momentum transfer. This originates technically from 
a cut-off of vacuum condensate terms in the short distance expansion where only the operators of 
lower dimension remain, and physically from neglect of the finite correlations between the quarks or 
the gluons in the physical vacuum. The case can improve by using the operator product expansion 
(OPE) around the light cone x 2 ~ instead of at the short-distance x ~ 0. The calculational 
framework established in such a way is termed light-cone QCD sum rule [3,4], which combines 
the well-developed description of exclusive process in terms of perturbative QCD factorization 
formula with traditional QCD sum rules. It has found wide applications [3-11] in the literature, 
especially applying to a study on heavy-to-light transitions, and the obtained predictions are 
encouraging. However, new nonperturbative inputs, the light-cone wavefunctions of light mesons 
classified by twists, are involved in light-cone sum rules. A question we must answer is whether 
these wavefunctions correctly reflect the underlying nonperturbative dynamics, in other word, 
whether they are phenomenologically better known. To date, only the components of leading twist 
have undergone a systematic investigation [6] for the distribution amplitudes of the nonsinglet, 
while the high twist wavefunctions, some of which play a role as important as the twist- 2 ones in 
most cases, require a further examination or refinement for improving the accuracy of numerical 
results. To prevent the resulting sum rules from suffering the possibly sizable contamination by 
the twist-3 components, an effective prescription || |TTJ has been suggested to estimate heavy-to- 
light form factors. The trick, to be specific, is to choose an adequate chiral current correlator, as 
allowed by the sum rules. As a result, not only the resulting new sum rules work effectively but also 
some of the twist-3 components, which make predominant contributions with respect to the other 
twist-3 wavefunctions in conventional light-cone QCD sum rule calculations, vanish automatically. 
From the theoretical viewpoint, this procedure offers a certain modification of the existing light- 
cone QCD sum rule approach to heavy-to-light transitions, making the nonperturbative effects 
well controlled. 

Also, among the important exclusive processes are the B decays with a photon emission. They 
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have in fact received a variety of model investigations. It is in [4] that the first application was 
presented of light-cone QCD sum rule technique to this type of processes, where the radiative 
decay £ — > P7 is discussed in detail. Later on, the same approach was applied to study B — > 
P7 [7,8], B — > /z/;7 [9] and B* — > B q ^ (q = u, dors) [10]. However, the problems with B* — > 
-B g 7 have been at issue. The yielded decay widths in different approaches disagree in numerical 
results and, to some extent, also in order of magnitude. For instance, light-cone QCD sum rules 
[10] predict T (B* u -> B u ^) = 0.63 KeV and T (B* d -> Btf) = 0.16 KeV while chiral perturbation 
theory [12] gives T (B* -> ^7) = O.UKeV and T (£3 -> 5 d 7) = 0.09 KeV. A definite conclusion 
depends strongly on our ability to calculate the relevant B*B q ^ coupling constant g_B*B 97 , which 
is intrinsically a nonperturbative quantity and is defined as 

7 (q, V) B* (p, e) I B q (p + q)) = gB*B ql e a(3p(T p a q fi r] p e a x . (1) 



where 77^ (e M ), g (p) are the polarization vector and the momentum of the photon(the B* meson) 
respectively. 

In the paper we would reexamine the radiative decay B* — > Bq'j by performing a calculation on 
the coupling Qb* B q ^ with an improved light-cone QCD sum rule approach, in which the appropriate 
chiral currents operators act as the interpolating fields. Our findings, are intriguing: The resulting 
sum rules may not only be of a clear physical meaning but also avoid the pollution by the nonlocal 
matrix element(7 (q)\ q (x) 7^75? (0) |0) starting with twist-3, which is known very poorly and thus 
may bring a large uncertainty to the sum rule results. 

This presentation is organized as follows. In the following section we propose our improved 
approach. Then a detailed light cone QCD calculation is preformed of the relevant correlator in 
Section 3. Section 4 is devoted to a numerical discussion on the sum rule results for (?_b*_b, 7 and 
then an estimate of T [b* — > B q ^ . We give a simple summary in the last Section. 

2. CHIRAL CURRENT CORRELATOR 

At the quark level, the radiative transition B* — > B q j is dominated by the magnetic moment 
interaction between the light quark and electromagnetic field due to small recoil of the decaying 
meson. This fact reveals, from the viewpoint of QCD, that the process in question belongs in the 
long distance type. It is well known that the current operators used as interpolating fields are by 
no means uniquely determined in QCD sum rules, and it is allowed to take different correlators 
providing that the validity of the resulting sum rule remains. This offers a possibility to explore 
the radiative decay B* — > B q ^ in a way consistent with the physical picture. In order to obtain 
such a desired sum rule result for gB*B q -y it is in order to adopt the following correlator, 



F, (V, (p + qf) = i J d'xe^ (7 (q, r,) \tJ^~ a (x) J s+p (0) 

= e^pfq^F (p 2 , (p + qf) + Vector Part, (2) 

where the chiral current operators J^~ A = qjfi (1 — 75) b and J^ +P = bi(l + 75) are chosen 
as the interpolating fields. For the invariant function F (p 2 , (p + q) 2 ), remarkably, the resulting 
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hadronic representation F H (p 2 , (p + q) 2 ^j depends not only on resonances of J p = 0~, 1~ but also 
on those of J p = + , 1 + . However, we can safely isolate the pole contribution of B q and B* , due 
to a sizable mass difference between the lowest negative and positive parity states, and express 
the high state contribution starting from the threshold parameters Sq and s' , which are set around 
the squared masses of the lowest 1 + and 0~ states respectively, in a form of dispersion integral. 
The result is 

F H (p 2 ( v + q) 2 = m ^ m ^ JbJb^b^ 

V A ' m b + m q (p2 _ m 2^ [(p + q f _ m 2^ 

I ( d SldS2 — p 4^4 — T , (3) 

J J (p 2 - Si) {p + q) —nip 



+ 



with the decay constants defined as 
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(B\bi l5 q\0)= ^ f Bq , (4) 
m b + m q 



[0\q ltl b\B;)=m Bq f B ^\ (5) 
After the Borel improvement P 2 — > M\ , (p + g) 2 — > M|, Eq. (3) is of the following form: 



'(6) 

3. LIGHT CONE OPE 

In this section, we will make use of the light cone OPE and give a detailed derivation of the 
theoretical expression for the correlator in Eq.(2). Schematically, the main contributions in the 
OPE come from the four Feynman diagrams shown in Fig. 1. All these diagrams can be calculated 
one by one for large space-like momentum regions p 2 and (p + q) 2 0, where the b quark is 
far off shell so that the OPE goes effectively. 

(1). SHORT DISTANCE CONTRIBUTION 

The short distance contribution, as will be shown, is less important from Fig.l (a) and (b) but 
its derivation turns out to be more complicated with respect to the long distance case. We have 
the UV-convergent loop integrals to deal with: 

1 



/ = / <fk 



\{p + k) 2 — m 2 ] [(k + p + q) 2 — m 2 ] (k 2 — m 2 ) 
k» 

[(p + k) 2 - m 2 } [{k +p~+q) 2 - m\\ (k 2 - m 2 ) ' 



p — [ d 4 k - (7) 

J \(p+k) 2 -m 2 ]\(k + p + q) 2 -m 2 ](k 2 -m 2 Y K ) 
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for Fig. 1(a) and the similar ones for Fig. 1(b). For the convenience of making the Borel transfor- 
mation, we use the forms of Feynman parametrization for them. After doing this, performing the 
loop integrals we have the invariant function F(p 2 , {p + q) 2 ), 

C 1 p 1 tTl T/Y. 

(8) 



i^V, (P + <lf ) = %[dx ^ dy _ m f-- r— 12-= 

Z7r^ jo j o m^a; + m z x — p z xx y — [p + q) xxy 



for Fig.l (a) and 



F {b \p 2 ,( P + q) 2 ) = 3 ^ f'dx fdy— 
2ir z Jo Jo mix 



m q x 



2tt 2 Jo Jo ° mix + m 2 x — p 2 xxy — (p + q) 2 xxy 
for Fig.l (b), with x = 1 — x and y — 1 — y. The total perturbative contribution can be written as 



(9) 



F (a+b) (p 2 ,(p + qf) 



3m q f 1 f 1 (Q q x + Q b x) x 



[ dx [ dy- 

Jo Jo j 



(10) 



27r 2 Jo Jo b m 2 x + m 2 x — p 2 xx y — (p + q) 2 xxy 
To express it in a form of dispersion integral we use the exponential form 

F {a+b \p 2 ,(p + q) 2 ) = ^ C dx C dy dax(Q q x + Q b x)e-<< w+m > +Q ' x ^ +Q ^ y ), (11) 

27T Z Jo Jo Jo 

where Q\ = —p 2 and Q\ = — (p + q) 2 . Following the method proposed in [13], we arrive immedi- 
ately at the perturbative spectral density p( a+fe ) (s±, s 2 ), 
3m„ r X2 dx 



P (a+b) (S!,S 2 ) 



with 



T^t J "=T ( m Q x + m b x ) 5 ( s i - s 2) © (si - (m b + m q ) 2 ^ (s 2 - (m b + 



m„ 



(12) 
(13) 



Xi = ^Si +m 2 b ~m 2 q - \js\-2 (m 2 b + mfj s x + (m 2 b - m 2 ) /2s u 

x 2 = ^si +m 2 b -m 2 q + )J s\ - 2 (ml + mf) Sl + (m 2 b - m 2 )^J j2s x . (14) 
Furthermore, completing the integral over x yields 

P {a+b) (si,s 2 ) = ^5(s 1 -s 2 )Q(s 1 -(m b + m q ) 2 )Q(s 2 -(m b + m q ) 2 ) 



x 



(Q6 -g,)A(M, B ) + Q,l„^4±|±Ai^| 



with A = ml/si, B = m 2 /s 2 and A (1, A, B) defined as 



A (1, A, B) = VI + A 2 + B 2 - 2A - 2B - 2AB. 
Putting everything together, we have the Borel improved result 



(15) 



(16) 



F [a+b) (Ml Ml 



^ m jdsQ(s-(m b + 



m n 



2\ -- 



e m 2 



x 



27r 2 Af 

in nun a />>, n i 1 ~ A + B + A (:L A D) 
{Qb - Q q ) A (1, A, B) + Q 9 In + ^ 



(17) 



+ B-X(1,A,B) 

with M 2 = MfMf/(Mf + M 2 ). It is explicitly shown that the short distance contribution is 
strongly suppressed from the perturbative hard quarks owing to smallness of the light quark mass 
m q , just as one may expect on the basis of a naive analysis. 



(2). LONG DISTANCE CONTRIBUTION 

Fig.l (c) and (d) denote the long distance effects induced by the light quark propagating in 
vacuum and interaction of the light quark in vacuum with the external electromagnetic field, 
respectively. In the light cone OPE, these soft quarks in vacuum are considered to be in a state of 
finite correlation, namely they are of a non- vanishing momentum [14], which stems from the high 
dimension operators neglected in conventional QCD sum rule calculations. We might describe the 
effect by introducing the concept of nonlocal quark condensate [14]. 

In the first place, we consider the contribution of Fig.l (c). In this case the nonlocal quark 
condensate (g(x)g(0))appears in the light cone OPE. To model its behavior near the light-cone 
x 2 = 0, a pragmatic strategy is to assume it to obey a distribution of the Gaussian type [14] 

(q(x)q(p)) = (qq)e-%, (18) 
with x\ = —x 2 . If comparing Eq.(18) with the short distance expansion 

2 

(q(x)q(0)) = (qq) ~ jl (qq) xl + ■ ■ ■ (19) 
with ml = 0.8 GeV 2 , the parameter p is found to be 5 GeV~ 2 . Introducing the parameters a, f3 
and 7 = a + (5 + p and then performing the cumbersome loop-integral, we obtain the following 
result for the contribution of the nonlocal quark condensate, 

F^{p 2 ,{p + q) 2 ) = -2p 3 Q b (qq) f° dad(5-^ P ^ Q --^ )m \ (20 ) 

Jo 7 d 

where P| = —p 2 ,Q 2 E = — (p + q) 2 - Furthermore, it can be converted into the Borel transformed 
form 

i Pm *( . M t +M l .. \ ( 1 I \ 
F" ) «,«3=-2«.r«) 5S « 4^«^J e (p- w - w ). (21) 

The O function restricts the Borel parameters M 2 and Mf to the range + < p. We have 
already checked that as p — > oo the above result returns precisely to what it should be in the case 
in which only the leading term remains in (19). 

A subtle approach has been suggested to deal with that type of dynamics in Fig.l (d). The 
idea is to introduce so called photonic light-cone wavefunctions to parametrize the long distance 
dynamics. This can be explained specifically as follows: On contracting the b quark operators, 
one has the matrix element T ~ iQ q r] ( ^ J dze iq ' z (0\JP m (z)O(x,0)\0), where J p em {z) is the electro- 
magnetic current operator of a light quark q with the charge Q q and 0(x, 0) expresses a certain 
nonlocal operator. If the underlying light quarks are hard we can perform the standard OPE. 
This is typically the case of Fig.l (b). For the soft light quarks we leave the light quark operators 
without contraction and treat J^ m {z) as the interpolating field of a photon state with the momen- 
tum q. In such a way we can define a matrix element of the nonlocal operator 0(x, 0) between 
the vacuum and the photon state as 

( 7 (?)|O(x,0)|0> = tQ4 x) j dz e*'(0 | J? m (z)O(x,0) | 0), (22) 
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and further parametrize it via a series of photonic light cone wavefunctions, which is equivalent 
to a summation over all the relevant condensate terms in the case of the short distance OPE, 
as in the applications of light cone QCD sum rules to B — > n transitions. With such a trick, a 
straightforward calculation of Fig.l (d) yields 



d A x J d 4 ke^~ k > x (7 (q) \q (x) ^ (1 + 75) q (0) | 0) 



Using the 7 algebraic relations 



mi 



k 7 



(23) 
(24) 



and considering the identities 



(7(?) I? W 755(0)|0> = 0, 
( 7 (g)|g(a;)g(0)|0) = 0, 



(25) 

(26) 
(27) 



we have 



*JV,(P + ?) 2 ) 



2i 



d 4 x J d 4 ke l{p - kyx 



k l 



rn 



2 _ k 2 [^7(9) \q(x)(T^q(0)\0) 



+ 



1 



W(7 (q) \q(x)a a/3 q(0)\0) 



{21 



An important observation is, as an appealing feature of our approach, that all the long distance 
dynamics entering the invariant function F^ d \p 2 , (p+q) 2 ) is encoded in the nonlocal matrix element 
(7 {q) \q ( x ) &apq (0)| 0), only one surviving after contracting the b quark operators. This answers 
very well to the actual physical picture in B* — > B q '~f since such a matrix element stands typically 
for an effect of the magnetic interaction. 

Light-cone expansion of (7(5) \q(x)a a /3q(0)\0) involves, to next-to-leading order in x 2 , the 
leading twist 2 wavefunction ip (u) and the two twist 4 distributions f\ (u) and / 2 (u) . The explicit 
form is 

(7 (?) \q (x) a aP q (0) I 0) = Q q (qq) Qf du\^p (it) F af3 (ux) + J dux 2 f x (u) F a/3 (ux) 



+ 



f duf 2 (it) 
Jo 



xpx^Fa^ (ux) — XaX^F^ (ux) — x 2 F a/3 (ux) J , (29) 



where (qq) denotes the quark condensate density, x indicates the magnetic susceptibility of the 
light quark and F a/3 is the external electromagnetic field tensor. Considering a simple situation 
q — > [4] is helpful to shed light on the physical meaning of Eq.(29). Substituting Eq.(29) into 
Eq.(28) yields 

ip(u) 



F {d \p\{p + qf) = 2Q q (qq)(x C du-^ 

\ Jo ml — 

- C du [/i (u) - f 2 (u) 
Jo 



(p + uq) 2 
4 



ml — (p + uqY 



+ 



^m\ 



m\ — (p + uqY 



.(30) 
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An emphasis we would put is that the contribution of the high Fock-state qqg, which should be 
taken into account to the present accuracy, has been neglected for it is expected to be negligibly 
small from most previous calculations. For later continuum subtraction, we have to convert Eq.(30) 
into a form of dispersion integral. We do that for the twist- 2 term and remain the twist-4 ones 
unaffected. The relevant spectral density is easy to evaluate by virtue of the technique in [13] once 
again. The calculation procedure is as follows. First of all, we apply a double Borel operator in 
p 2 and (p + q) 2 to the twist- 2 term, obtaining, 



S(MM)B(Ml,«)jf'* *j 



M l + M 2~2 



(p + uq) 2 M\ + Mf 



M\ 



Ml + Mi 



(31) 



The symmetry of the correlator in Eq.(2) allows us to set Ml = Mf so that the wavefunction 
x can take its value at the symmetric point /1q = 1/2. Considering it and making a 



replacement Mf — > Mf — > ^- in Eq.(31), we get 



B(-,Qi)B[-,Q 



0-2 



o m 2 — (p + uq)' 



_^_e-* +CT2 V(l/2) 

0\ + 0" 2 



= / (0-1,0-2) • 



(32) 



Finally, we take the function -^-J (oi,o 2 ) = / [pxp-i) and perform one more Borel transformation 
in the variables U\ and a 2 . The relevant spectral density reads 



(tw2) 



(si,s 2 ) 



sis 2 



B[-,a 1 )B(-,a 2 )f(a 1 a 2 ) 



Sl 



\s 2 



= S( Sl - s 2 ) y?(l/2)0 ( Sl - m 2 ) (s 2 - m 2 b ) . 



(33) 



At present, it suffices to represent the twist-2 contribution in a form of dispersion relation. We 
derive the following expression for F <yd \p 2 , (p + g) 2 ), 



F^(p 2 , (p + q) 2 ) = 2Q q (qq) [ X I d Sl ds 2 



5( Sl -s 2 )ip(l/2) 







du [/j (u) - f 2 (u)} 



(p - si) \{p + q) - s 2 
4 



8m 2 



m 2 - (p + uq) 2 ml — (p + uq) 



.(34) 



After the Borel transformation it becomes 

-=(0 f -. ,0 . -o\ i 



(m 2 ,m 2 2 ) 



Ml Ml 



2Q q (qq) V (1/2) x / , dse~^ - 8Q q (qq) 



mi 



x [/ l( l/2)-/ 2 (l/2)](l + ^ e" 



(35) 



Now we conclude this section with writing down the light-cone OPE result for the invariant 
function F (p 2 , (p + q) 2 ^j in the Borel transformed form, 



p(OPE) 



(Ml, M, 



(Ml, Ml) + F (c) (Ml, Ml) + F (d) (Ml, M 



(36) 
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4. SUM RULES AND NUMERICAL DISCUSSIONS 



Matching Eq.(36) with the corresponding hadronic representation (6) and using the quark- 
hadron duality ansatz, in which the hadronic spectral density p H (si,S2) is usually assumed to 
coincide with that derived in the light-cone OPE, we arrive at the sum rule for the product 

fB q fB*gB*B ql , 

m b + m q B _i B „ _ ^ _jo \ 2 



fB q fB*9B*B ql = J __ q e ^ <| 2Q q (qq) X y? (1/2) - j M - 2Q b (gg) e ^ 



m%m B * 



M- 



/ x 1 — m?/s + m^/s + A (1, mils, mils ) 

(Q6 - Qq) A 1, mg/ a , mils) + Q q In ' J 

1 — mys + my s — A 1 1, mys, my s 



x 

1 - m^/s + m^/s - A (1, m£/s, ra^/sj 

(37) 

At this point, we have a minor comment. Nonperturbative QCD dynamics dominates the sum 
rule; especially the magnetic interaction term makes the leading contribution, while the hard quark 
contribution modifies merely the sum rule at the level of m q , a negligibly small effect. The nonlocal 
matrix element (7 (g) \q (x) j^q (0)| 0), which starts with twist-3, is effectively eliminated, making 
the sum rule free of the contamination by its uncertainty. In addition, what we should emphasize 
is that our subtraction procedure is different from that in Ref. [10], in which continuum subtraction 

2 2 
rn b m b s 

e~M^ — > e~M? — e^M 1 is made for the twist-3 and -4 parts. 

The actual calculation needs a dynamical input. The nonperturbative parameters concern- 
ing the light quarks include the light-cone wavefunctions <p(l/2), fi (1/2) and / 2 (l/2), quark 
condensate (qq), magnetic susceptibility x an d light quark mass m q , Among these, the leading 
twist-2 wavefunction ip {1/2) embodies the underlying dynamics and hence demands a precision 
determination. Fortunately, it follows from a well-found calculation that ip (u) is dependent quite 
weakly on the scale and deviates little from its asymptotic form [4] , with a well-tried universality. 
Therefore, we can use 

tp (u) = 6u (1 - u) , (38) 

to high accuracy. According to the analysis given in [7], the set of the twist-4 wavefunctions can 
be specified as 

/i(«) = ~(l-«)(3-«), (39) 

h{u) = -\{l-uf. (40) 

Since we work in the case that the QCD radiative correlations are neglected, the appropriate 
scale, at which the quark condensate (qq) and magnetic susceptibility \ t a ke values, should 
be set by the typical virtuality of the underlying b quark fib ~ V m B q ~ m h as m Ref. [8]. 
We take the results used widely (uu) (/i = 1 GeV) = (dd) (// = 1 GeV) = - (0.24 GeVf , 
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(ss) (/i = 1 GeV) = 0.8 (uu) (/i = 1 GeV), and x(» = 1 GeV) = —4.4 Ge\/- 2 [15], from which the 
use of the Renormalization Group Equation (RGE) gets (uu) (/it) = (^dctj (yU^) = —0.018 GeV 3 , 
(ss) (fib) = —0.014 GeV 3 and x (t*b) — —3.4 GeV~ 2 . Furthermore, the light-quark masses are set 
as m u = irid ~ 0, m s = 0.115 GeV. As far as the B channel parameters are concerned, we use 
nib = 4.8 GeV for the b quark mass and tub a = 5.279 GeV, rriB* = 5.325 GeV, rriB = 5.369 GeV 

u,a u ^ * 

and ms* = 5.416 GeV for the set of the B meson masses. The relevant decay constants, for con- 
sistency, require a recalculation in the two point sum rules in which a chiral correlator should 
be chosen in a proper way. Moreover, we would do that at the leading order in a s (Q 2 ) since 
in present case the sum rule for f B q f B* q 9B* q B q ~i works at the same accuracy. In the light of the 
prescription in [11], the best fits lead to f Bud = 120 MeV, f B * d = 141 MeV with s = 32 GeV 2 , 
and f Bs = 147 MeV, f B * = 160 MeV with s = 34 GeV 2 . All the above parameters will be used 
as the central values in the following numerical estimates. 

To proceed, one has to find the reasonable ranges of M 2 , from which the desired sum rules can 
be read off. This is a critical step towards deriving a reliable prediction. For this end, we use the 
standard procedure where the high state contributions are limited to a level below 30% while the 
twist-4 ones to an order less than 10%; simultaneously it is required that the resulting sum rules 
be considerably stable. In the first place, we focus on the derivation of the sum rule for g B *B u -y In 
this case, using the above criterion the fiducial values of M 2 are found to lie between 8 — 14 GeV 2 , 
with so = 32 GeV 2 . In this range, M 2 dependence of the product /s tl /s;5 , B*B ll7 , as shown in Fig. 
2, is quite weak. From the "window" we have the sum rule results 

f B jB Z 9B Z B^ = 10.12 ± 0.17 MeV, g B * Bul = 0.59 ± 0.01 GeV' 1 , (41) 

the quoted errors being due to the change of M 2 . Furthermore, it is necessary to look into the 
uncertainties arising from the input parameters. Let Sq change between 31 — 33 GeV 2 while 
the other parameters keep fixed, the resulting variation of the product fB u fB*gB*B u -y is roughly 
±8%. Utilizing this result and incorporating the corresponding changes of the decay constants, 
the uncertainty in the coupling gB*B u -y is estimated at the level of ±5%. To look at the impact 
of the uncertainty in mj,, we concentrate on the sum rule for gB*B u -y and consider a correlated 
variation effect, which appears as changes from 4.6 GeV to 5.0 GeV , in such a way in which 
we substitute the analytical forms of the sum rules for the decay constants into Eq.(37), observe 
the variation of gB^B u -y with by requiring that all the sum rules take only their values from the 
best fits in both sq and M 2 . A detailed analysis shows that in this case the resulting influence 
amounts to ±6%. We add also a numerical uncertainty of 30% to all the condensate parameters 
and see the effects. This results in an error of ±8% in the sum rule for gB*B u -y- Finally, the 
resulting total uncertainty in gB*B^ can be fixed at about 21%, by adding linearly up all the 
considered errors. Also, the same procedure is applied to the case of B* d B,pi; however, we would 
like to consider an additional source of uncertainty in the B* S B S ^ case, allowing m s to change in 
the region of 115 ± 25MeV. The resulting two sum rules turn out to have the Borel intervals and 
uncertainties different slightly from the corresponding those in the sum rule for fB u fB*gB*B u 'y- We 
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don't give details any more for brevity. The final results read: 

fB d fB*gB*B dl = -5.31 ± 0.15 MeV, g B * Bdl = -0.31 ± 0.01 GeV' 1 , (42) 
with s = 32 GeV 2 and M 2 = 7 — 14 GeV 2 , and 

f B jB* s gB t B sl = -7.46 ± 0.23 MeV, g B . B „ = -0.32 ± 0.01 GeV' 1 , (43) 

with s = 34 GeV 2 and M 2 = 8 — 15 GeV 2 . The total uncertainties in the two coupling constants 
are respectively about 20% and 23%. The sum rule stability is illustrated in Fig.2 too, for the 
products f B J B *g B * Bin and fs a fB* 

Having the sum rule results for g B * Bql at hand, we can calculate the decay widths by means 
of the formula 

^ g qJJ 96n [ m B * J 

The results are 

r(5* -> B ul ) = 0.89 ± 0.34 KeV, 

Y(B* d -> B d -i) = 0.25 ± 0.10 KeV, 

T(B* -> B sl ) = 0.28 ± 0.13 KeV 

It is obvious that the sum rule results (45) and (46), within the errors, are in agreement with 
the estimates from the standard light-cone QCD sum rules [10]: T (B* — > B u ^y) = 0.63 KeV and 
T (B* d — > B(Pi) = 0.16 KeV . If we make a comparison of our sum rule predictions and chiral 
perturbation theory ones [12], on the other hand, it is demonstrated that there is a considerable 
numerical difference in the obtained decay widths in the case of B* — > B u ^; but the obtained 
predictions are comparable with each other within errors in the case of B d — > B<ff. Introducing 
the ratio R — T (B* — > B s ^) /T (B' d — > B d ^y) to see SU (3) breaking effect in the radiative decays, 
we find R = 0.43 — 2.72, with a large uncertainty. However, the resulting lower limit of the ratio 
is close to R — 1/3 0.33 obtained in chiral perturbation theory [12]. 

5. SUMMARY 

To sum up, in this paper we present an improved light-cone QCD sum rule approach to the 
radiative decay B* — > Bcpi whose model estimates are at issue, an adequate chiral current correlator 
being used to get a sum rule result consistent with the underlying physics. In this framework, 
a detailed derivation is given of the sum rules for the relevant coupling constant g B * Bql , which 
parametrizes all the long distance dynamics. The resulting sum rules are of the two main features: 
(i). It is the magnetic interaction which makes the leading contribution to the resulting sum rules; 
on the contrary, the perturbative hard quarks modify the sum rules only at the level of m q . This 
typically is the principal dynamical feature the radiative decay B* — > B^ possesses, (ii). The 
nonlocal matrix element (7(g) \q (x) 7^75? (0)| 0), which is equally important but poorly known 
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(44) 

(45) 
(46) 
(47) 



in comparison with the leading matrix element (7 (q) \q(x) u a pq (0)| 0) [10], disappears in our 
sum rules, which controls effectively the pollution due to the long distance effects. In addition, 
we apply the light cone OPE instead of the short distance expansion to the quark condensate 
contribution, minimizing the uncertainties in the sum rules to the accuracy in consideration and 
guaranteeing a self-consistent sum rule calculation. It is predicted that V (B* — > 5«7) = 0.89 ± 
0.34 KeV, T (B* d -> B dl ) = 0.25 ± 0.10 KeV and T (B* s -> B s -f) = 0.28 ± 0.13 KeV All these 
channels are predicted to be of the decay widths at the level of 10~ 7 GeV and therefore are 
accessible at the current running B factories and the future LHC. Also, a comparison is made 
with other model predictions. There is an approximate agreement between our results and those 
from the standard light-cone sum rules [10], for V (B* — > B u ^) and V (B d — > B^). We find also 
that the resulting sum rule predictions deviate more or less from those of chiral perturbation 
theory [12], depending on different cases. Striving for perfection of phenomenological models 
and enhancing ones' ability in controlling nonperturbative dynamics, we believe the procedure 
presented here to be a valuable attempt in the direction. 

Z. H. Li thanks Prof. K. T. Chao for helpful discussions and comments. This work is in part 
supported by the National Science Foundation of China. 
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(a) 



(b) 




Fig.l. Diagrams making the predominate contributions to the correlator (2). (a) and (b) depict 
perturbative contributions, (c) nonlocal quark condensate correction and (d) emission of photon at 
long distances parametrized by photon wavef unctions. Solid lines denote quarks, crosses vacuum 
quark condensates and wavy line photon. 
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Fig.2. The light-cone QCD sum rules for f Bq f Bq g B * BqJ , with s = 32 GeV 2 for f Bu f B * u g B * uBul 
(dotted) and f B jB*gB*B dl (dashed), and s = 34 GeV 2 for f Bs fB* s gBiB sl (solid). It should be 
understood that we have set f Bd f B * d g B * dBdl and f Ba f B *g B * Bsl positive, for a comparison. 
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